Fractional charge excitations in fermionic ladders. 
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The system of interacting spinless fermions hopping on a two-leg ladder in the presence of an 
external magnetic field is shown to possess a long range order: the bond density wave or the 
staggered flux phase. In both cases the elementary excitations are Z2 kinks and carry one half the 
charge of an electron. 
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It has long been established both theoretically Q, Q 
and experimentally 3] that quantum numbers of elemen- 
tary excitations in interacting many-body systems are 
not necessarily limited to the values characterizing free 
particles. While excitations carrying, for example, an 
electric charge double the free electron charge e can be 
visualized as pairs of the original particles |4j , the frac- 
tionally charged excitations may appear counterintuitive 
Li LJ LJ] By far the most celebrated example of such exci- 
tations are the fractionally charged quasiparticles in the 
quantum Hall state At about the same time the no- 
tion of fractional charge appeared in the Peierls model ap- 
plied to quasi-one-dimensional conducting polymers |2|. 
More recently [fjj, fractionally charged excitations were 
predicted in the context of the extended Hubbard model. 

In this paper we argue that charge fractionalization 
is a generic phenomenon in the sense that it appears al- 
ready in the simplest model of spinless fermions in a mag- 
netic field. Restricting the fermions to hop on a two-leg 
ladder0,0, we show that just nearest-neighbor interac- 
tion (if it is strong enough) leads to long range order 
(LRO) with a doubly degenerate ground state. The el- 
ementary excitations are then quantum domain walls or 
Z2 kinks that carry the charge e/2. 

The applied field plays a crucial and somewhat surpris- 
ing role in the problem. In the absence of the field the 
only possible LRO is the bond density wave (BDW) sim- 
ilar to that of Ref. 0. However in the presence of the field 
different types of LRO are possible. The most spectacu- 
lar manifestation of the field is the staggered flux phase 
[9L flo| illustrated in Fig.^ Note that the stag gere d order 
parameter here is induced by a uniform field . 

To build some intuition about the ordered phases con- 
sider the strong coupling limit of the model where the 
transverse interaction V± is the strongest, so that no rung 
of the ladder can be doubly occupied. For the quarter- 
filled ladder the starting point is a state where there is 
one electron per two plaquettes. The hopping terms tend 
to delocalize the electrons along the links. If the longitu- 
dinal interaction is repulsive Vj| > 0, then the most favor- 
able configuration is that in which the electrons, avoiding 
to reside on neighboring sites of the same chain, get delo- 
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FIG. 1: The staggered flux phase in the two-leg ladder pene- 
trated by a magnetic field. The arrows inside the plaquettes 
indicate the direction of the applied uniform field (directed 
inside the plane of the Figure). The arrows along the links 
indicate directions of the staggered currents. The central pla- 
quette is the domain wall between the two degenerate ground 
states. This excitation carries the charge e/2. 



calized on rungs thus leading to a commensurate BDW. 
If, on the other hand, V\\ < then the system has a 
natural tendency to phase separate which is opposed by 
the flux. If the flux is strong enough the electrons are 
placed on every other plaquette and are completely de- 
localized around them (thus avoiding double occupancy 
of the rungs) producing currents circulating in the same 
direction (determined by the field). This can be viewed 
as a staggered current superimposed upon the average, 
uniform (persistent) current. 

In this paper we show how the above cartoon picture 
arises starting from the weak coupling limit. Our strat- 
egy is the following: (i) we diagonalize the Hamiltonian 
in the absence of interaction; (ii) linearize the exact spec- 
trum near the Fermi points and bosonize the model; (iii) 
finally, we solve the bosonic problem and find the ground 
state of the interacting system. 

The single-particle spectrum of the problem consists of 
two one-dimensional bands and is governed by the mag- 
netic flux through a plaquette in units of the flux quan- 
tum / as well as the ratio of the transverse and longi- 
tudinal hopping amplitudes r. It is possible to partially 
fill only one band and have only two Fermi points. In 
this case the low energy physics of the problem is similar 
to that of the XXZ spin chain ^^l- If the interaction is 
weak, then the system is in the Luttinger liquid regime 
(i.e. similar to the gapless phase of the XXZ chain). If, 
on the other hand the two Fermi points are commensu- 
rate with the lattice and the interaction is sufficiently 
strong, then the umklapp scattering becomes relevant in 
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the renormalization group sense and opens a gap. The 
ground state of this gapful phase possesses LRO and is 
doubly degenerate. Now, due to the applied field, the 
physics of the ladder is richer than that of the single 
chain. Depending on the sign of the umklapp term dif- 
ferent types of LRO are possible: BDW in the repulsive 
case, or the stag gered flux phase [j| (or an orbital anti- 
ferromagnet |1JJ]) in the case of attraction. 

Having sketched our line of reasoning we now describe 
our calculations. In this paper we focus on the particular 
parameter regimes where the system exhibits the frac- 
tional charge. Mathematical details and a description of 
the full phase diagram of the system as well as a quan- 
titative description of the strong coupling limit will be 
discussed elsewhere [lflj . 

We start with the tight-binding Hamiltonian 

n = -^2[\ (t\\(yi)4( x n)ci(x n +l) + h.c.) 



i=l,2 



-t±c{(x n )c 2 (x Tl 



h.c. 
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(1) 



where Cj(x n ) is the electron annihilation operator on the 
chain j at the site x n , ij_ and t\\ are the transverse and 
the longitudinal hopping amplitudes. The magnetic field 
B is introduced by means of the Peierls substitution . 
Choosing the Landau gauge |l5| with the vector potential 
A = B{— y,0,0) and defining the y coordinates of the 
chains as 2/1(2) = ±6/2 we write the longitudinal hopping 
amplitude as 



*l|(f)=*| 



(2) 



where / = Bab/<po, </>o being the flux quantum and a the 
lattice spacing along the chains. The control parameter r 
is defined as r = t±/t^; tff' plays the role of the band- 
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width and hereafter will be set to unity. Expressed in 
terms of the flux the model is explicitly gauge invariant. 




FIG. 2: The spectrum in the region r > 1, r cos7r/ > sin irf. 

The physical quantities we shall study in this paper are 
the electrical current and particle density. In terms of the 
lattice operators Cj(x n ) the in-chain current is given by 



Ji 



(2) 



"2 cJ (2) (x„)ci (2) (x„ + i)e =F7r,/ 
while the transverse current is defined as 
Jj_ = -it c\(x n )c 2 (x n ) - h.c 



h.c. 



(3) 
(4) 



Similarly, the bond density on transverse links is 

p± = c\(x n )c 2 (x n ) + h.c. (5) 

In the absence of interaction the Hamiltonian Eq. JI} can 
be diagonalized with the help of the linear transformation 

ci(q) = u q a q + Vgf3 q ; c 2 (q) = v q a g - u g f3 q ; (6) 

where the "coherence factors" are (the signs are explicit 
in Eq.© so that the coherence factors are positive) 



u2 q( V q) 



it 



sin q sin irf 



V sin 2 q sin 2 irf + t 2 



(J) 



The exact spectrum of the system consists of two one- 
dimensional bands 



cos k cos irf T- y sin k sin irf + t 2 



(8) 



where k in the momentum along the chains. Since the 
transformation / — ► 1 — / preserves the spectrum we 
only need to consider the fluxes such that < / < 1/2. 
If the transverse hopping amplitude is large enough so 
that t > cos irf the spectrum possesses a band gap and 
resembles a band insulator, see Fig. [21 If the flux is not 
too small, sin 2 irf > r cos irf, the bands acquire a double- 
well shape shown in Fig. [3] 




FIG. 3: The spectrum in the region (sin7r/) 2 /r > cos7r/ > r. 
The horizontal line depicts the chemical potential such that 
there are only two Fermi points. 

Now we focus on the situation where the lower band 
is only partially filled (the upper band is empty) and 
there are only two Fermi points. The ground state in 
the absence of interaction is characterized by the persis- 
tent current flowing along the chains. Indeed, the ground 
state expectation value of the transverse current Eq. (gj 
vanishes (J±) = 0, while the in-chain current Eq. |J3J is 
non-zero as long as the flux is applied 



(Ji) 



sin7r/ 



E 

q=-k F 



sin q cos irf 
\J sin 2 q sin 2 irf + r 2 



cos q 



.(9) 



The current flows along the two chains in opposite di- 
rections (Ji) = — (J 2 ), so that the net current is equal 
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to zero as it should. The Fermi momentum in Eq. © 
depends on the applied field as 

cosfcp = — jl cos ir f — \J t 2 + (1 — jl 2 ) sin 2 nf, (10) 

where jl is the chemical potential in units of tu* . 

It is noteworthy that all (filled) states contribute to 
the persistent current Eq. Q. The persistent current in 
this problem is not an infrared phenomenon and thus can 
not be accounted for by an effective low energy theory. 
The low energy effects we discuss below are thus taking 
place on the background of this single-particle persistent 
current. We will not discuss further corrections to Eq. 
in this paper. 

To formulate an effective low energy theory we now 
linearize the exact spectrum in the vicinity of the Fermi 
points Eq. 1)10(1 . assuming that that the k = features 
of the exact spectrum are far away from the Fermi level, 
the condition that has to be checked later. Since only one 
band is filled, this step is fairly standard and is similar to 
the case of a single chain|l2j. Separating fast and slow 
variables, we introduce the left- and right-movers for the 
lower band and find the standard kinetic term 



Ha 



(a R (k)a R (k) - al(k)a L (k)) , (11) 



where the Fermi velocity is 



Vp = sin k p 



\f{l — jl 2 ) sin 2 n f + T ' 
\J sin 2 fcp sin 2 irf + r 2 



(12) 



To obtain the final form of the effective theory we need 
to express the physical observables in terms of the low 
energy fields (in doing so we approximate the coherence 
factors by their values at the Fermi points uq and Vq). 
The smooth part of the particle density is then given by 



Pi = u J R + v J L ; p 2 = v J R + u J L , 



(13) 



where Jr(l) ='■ a R iL) a R(L) '■■ Mostly we shall be in- 
terested in staggered operators with expectation values 
vanishing in the absence of interaction. The operators of 
interest are the staggered currents (defined on links) 



(—1)" / 
3i{i){x n ,x n+ i) = — - — u Q VQ\a\(x n )a R {x n+ i) 



+i)e 



-iirf 



-a\{x n+ i)a R (x n )e™ f - h.c 



)■ 



(14a) 



3± 



(x n ) = it±(-l) n (yl - it 2 ) (a\(x n )a R (x n ) - h.c.^j , 

(14b) 

and the staggered bond density 



0) 



(x n ) = (-1)™ (a[(x n )a R (x n ) + 



(15) 



Current conservation requires \ (j±) \ — 2|(ji( 2 ))|. 

We now turn to the discussion of interaction effects. 
For simplicity we consider the nearest-neighbor density- 
density interaction 

Ti-mt = ^ \ V \\ \J l 'Li x n)ni(x n+1 ) + n 2 (x n )n 2 (x n+1 )^ 



+Vj_n 1 (x n )n 2 (x n ) 



(16) 



with rij —: CjCj : (colons indicate normal ordering). 

In terms of left- and right-movers (neglecting trivial 
renormalizations of the Fermi velocity) we rewrite the 
interaction Eq. (|16|) in position space as 

Hi n t » ^2 [31 : a R {xi)ot R {xi)a\{xi)a L (xi) : (17) 

i 

+ {g 2 : a R (xi)a R (x i+ i)a L (xi)a L (x l+ i) : +/i.c.J ; 

where the new interaction constants are [l6| 

.91 = 2V\\u 2 vl (1 - cos2/c F ) + V x {u 2 - v 2 f; (18a) 



52 = -2Vj,e 



(18b) 



The first term in Eq. (|17fl is the density-density interac- 
tion that can be taken into account exactly by means of 
bosonization. The second term is the umklapp interac- 
tion that does not conserve momentum and thus is only 
relevant when kF is commensurate with the lattice (i.e. 
for the quarter-filled ladder). 

In the commensurate case the Fermi mo mentum kp = 
ir/2 and the chemical potential jl = —y/ sin 2 nf + r 2 . 
This is consistent with the requirement of having only 
two Fermi points when 2r > | cos27r/|/cos7r/. In this 
case the combinations of coherence factors in Eqs. 1(17(1 
and 1(14(1 simplify to 

u l ~ v o = ( sin7r /)/A ; u v = T/fj,. (19) 

Now we use the standard bosonization procedure and 
express the interaction Hamiltonian Eq. 1(171) in terms 
of a bosonic field 4> and the conjugate momentum II. 
Combining it with the kinetic term we arrive to the usual 
form of the Luttinger liquid Hamiltonian (here a is the 
bosonic cut-off) 



n, U F 



ku 2 + i (d^y 



92 



2n 2 a 2 



■ cos • 



(20) 



where the Luttinger liquid parameter K is determined 
by the interaction constant Eq. I(18a(l as K — [(1 — 
g 1 /2nv F )/{l+g 1 /2nvp)] 1 / 2 . 

The behavior of the effective theory Eq. 1(20(1 is well 
known 12]. At the critical value of the Luttinger liquid 



4 



parameter K = 1/2 [lj| the s yste m exhibits a Berezinski- 
Kosterlitz-Thouless (BKT) |l8| transition to a gapful 
phase with broken Z 2 symmetry and LRO. The type of 
ordering depends on the sign of the coupling constant <?2 : 
(i) If the umklapp term is repulsive 52 > then 
the local minima of the potential arc achieved at <f> = 
\/tt/4, n = 0, ±1, In this case the stag- 



gered bond density Eq. ljl"5|) ~ smy/4Tr<fi acquires a 
non-zero expectation value and LRO is of the BDW-type. 
This ordering resembles the charge density wave in the 
extended Hubbard model Q at 1/4 filling. 

(ii) If, on the other hand, 92 < [l9(, then the local 
minima are at <f> = and now there exists the 

staggered current Eq. (|14fl j± ~ cos \/4ir(f) leading to the 
staggered flux phase illustrated in Fig. ^ Clearly, the 
currents are non-zero only in the presence of the applied 
field, as jj_ oc sin7r/. This is a surprise: we find that the 
uniform field causes the staggered current. 

In both cases the ground state is doubly degenerate 
and the elementary excitation is a "quantum domain 
wall" or, in terms of the effective theory Eq. (|2U|) . a Z2 
kink interpolating between the two neighboring minima. 
The charge carried by the kink is related to the distance 
A<fr between the minima [as follows from Eq. i|13|) ] 



Q 



Acj) 



(21) 



Since both types of ordering originate from the same in- 
teraction term in Eq. I|20() , the distance between the min- 
ima is the same for both cases, A<j) = \/7r/4, which brings 
us to the conclusion that the kinks carry fractional charge 



Q 



e 

± 2- 



(22) 



This statement is the main result of the present pa- 
per. Three features of the model were necessary to ob- 
tain this result: (i) the electron density is bosonized 
as p = dxip/y/n; (ii) there were two Fermi points with 
2kp = 7r; (iii) the interaction was strong enough to open 
the umklapp gap and produce the LRO. 

To summarize, we have considered the simplest possi- 
ble model of interacting electrons subject to an external 
magnetic field restricting the motion of electrons to hop- 
ping on a two-leg ladder. Despite its apparent simplicity, 
the model exhibits phase transitions to non-trivial LRO 
phases, the BDW phase in the case Vji > and the stag- 
gered flux phase for V\\ < 0. In both cases (due 
to their topological nature) elementary excitations (i.e. 
kinks) carry fractional charge. As such, charge fraction- 
alization appears to be a generic phenomenon mostly in- 
sensitive to the details of the model. As long as the basic 
features of the model (i.e. the transverse hopping and the 
flux) are in place any interaction that contains umklapp 
processes |l(j will lead to charge fractionalization. 



The applied magnetic field played a crucial role in our 
considerations. It led to a non-trivial deformation of the 
single-particle spectrum (Fig. OJl previously unaccounted 
for and to the appearance of the staggered flux phase. 
The Fermi momentum became strongly dependent on the 
flux. For a situation close to half-filling this implies the 
existence of four essentially different Fermi points. This 
and other features of the model that were not discussed in 
the present paper, as well as further generalizations of the 
model to larger number of chains and spinful electrons is 
the subject of a future publication |l3j . 

We are grateful to N. Andrei, T. Giamarchi, and A.M. 
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